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Abstract. We present CO(1-0) and CO(2-l) line emission maps of the barred spiral active galaxy HE 1029-1831 
(z=0.0403) obtained with the IRAM Plateau de Bure Interferometer (PdBI) and in part by the Berkeley-Illinois- 
Maryland Association (BIMA) observatory. Fhe CO emission is well associated with the optical bar and extended 
along it. Fhe FWHM of the CO emission is estimated to be ~(6±2) kpc. Fhe CO emission shows a strong velocity 
gradient along the minor axis of the bar (PA=90°). Fhe molecular gas mass is estimated to be ~1.2xl0 10 .M© 
which indicates a very gas rich host galaxy. Most of the molecular gas appears to be subthermally excited and 
cold but we also find weak evidence for a warmer and/or denser gas component at the southern part of the bar 
emission, about ~4 kpc from the galactic nucleus. 
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1. Introduction 

HE 1029-1831 is part of a nearby QSO sample (Bertram 
et al., in prep.) consisting of ~100 sources from the 
Hamburg/ESO survey for bright QSOs. The sole selection 
criterion of these sources is their low redshift of z < 0.06 
providing the capability of observations on the smallest 
observationally available angular scales and the accessibil- 
ity of several important diagnostic lines in the NIR (e.g., 
the CO (2-0) rotation vibrational band head absorption 
line) . We aim at determining the distribution and dynam- 
ics of molecular gas in the inner 1 kpc of these objects. 
This will add to a systematic study of the different mech- 
anisms of gas fueling into active galactic nuclei (AGN). 
Together with the NUGA PdBI survey containing Seyfert 
and LINER galaxies with z<0.013 (Garcia-Burillo et al. 
2003a,b,2005; Combes et al. 2004; Krips et al., 2005) this 
study will cover a large part of the activity sequence of 
AGN. Especially, the role of nuclear embedded bars, rings, 
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spirals as well as mini-spirals, warps and micro-warps in 
the fueling process is an important issue to investigate 
possible links between different nuclear disk morphologies 
and activity types. 

Thirty objects of this sample were observed in CO 
with the SEST telescope and ten with the BIMA array 
(Bertram et al., in prep.) resulting in a low detection 
rate of~15% which might be a result of a lack of instru- 
mental sensitivity and uncertainties in the given redshifts; 
these two observing campaigns were initially conducted 
to identify the objects with the strongest CO emission. 
As a consequence of the low detection rate, we have car- 
ried out deeper integration of an additional 26 sources 
with the IRAM 30m telescope leading to 19 further detec- 
tions (Betram et al. 2006, Bertram et al. in prep.). The 
BIMA observations resulted in a clear detection of the 
source HE 1029-1831 and in a tentative one of HE 1136- 
2304. Both have thus been observed as a follow-up with 
the IRAM PdBI array but only the first one could be de- 
tected and mapped while the tentative BIMA detection 
of HE 1136-2304 could not be confirmed. HE 1136-2304 
most likely lacks an exact redshift, too. 

HE 1029-1831 contains (m=2) spiral arms and a 
prominent bar (PA~30°) in the optical (H-band; Fischer 
et al. 2006; see Fig. [1]). Its redshift was recently deter- 
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Fig. 1. H-band image of HE1029-1831 (from Fischer 
et al. 2006) superimposed with contours of the inte- 
grated CO(2-l) line emission. Contours are from 3a =1.4 
Jy/beam km s _1 to 13 Jy/beam km s _1 in steps of 3cr. 



mined to be z = 0.0403 ± 0.0001 (Kaldare et al. 2003). 
The classifications of its active nucleus range from an HII 
region and extreme starburst to an AGN (Kewley et al. 
2001). However, recent NIR observations have given evi- 
dence that HE 1029-1831 may be classified as a luminous 
infrared galaxy with a weak narrow line Seyfert nucleus of 
type 1 (Fischer et al., 2006). 



2. Observations 

2.1. BIMA 

HE 1029-1831 is one of 10 members of the sample that 
were scanned for CO(l— 0) emission with BIMA at Hat 
Creek, CA, USA between October 2000 and April 2001. 
The interferometer, consisting of 10 antennae with a di- 
ameter of 6.1m each, was used in C configuration. The 
correlator was configured to provide 8 spectral windows 
in the upper sideband with 32 channels each and cov- 
ered an overall bandwidth of 800 MHz. The observations 
were scheduled over a range of several months with varying 
weather conditions, total integration time and uv-coverage 
for the different sources, hence resulting in different spatial 
resolutions and detection limits. In the case of HE 1029- 
1831, the BIMA survey led to a detection of bright CO line 
emission (Figj2|). These data were reduced using 3C273 as 
bandpass and 1048-191 as gain calibrator. A medium spa- 
tial resolution of 13.6" x 5.8" was achieved in this antenna 
array configuration. 
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Fig. 2. Integrated CO(1-0) emission as observed with 
BIMA (grey scale + dotted grey contours). The PdBI 
CO(1-0) (solid black contours) are superimposed on the 
image. The contour lines for the BIMA data represent mul- 
tiples of la levels, starting at 3c, where ler=3.0 Jy/beam 
km s . The contours for the PdBI CO(1-0) emission run 
from 3a —0.8 Jy/beam km s _1 to 7.3 Jy/beam km s _1 in 
steps of 3a. 

2.2. PdBI 

As BIMA follow-up, observations of the CO(1-0) and 
CO(2-l) line emission in HE 1029-1831 were conducted 
in February 2002 with 6 antennae in B and C configura- 
tions of the IRAM PdBI and again in March 2003 with 6 
antennae in A configuration. The bandpass calibrator was 
set to 3C273 and amplitude and phase calibrators were 
1048-191 and 1055+018. The latter were observed every 
twenty minutes. The 3.5 mm receiver was tuned to the 
redshifted 12 CO (1-0) line and the 1.2 mm receiver to the 
redshifted 12 CO(2-l) line. The amount of water vapor is 
calculated to be <6 mm for both sets of observations thus 
providing reasonable weather conditions and system tem- 
peratures for a low declination target. At both frequencies 
we used a total bandwidth of 580 MHz with a frequency 
resolution of 1.25 MHz. The phase tracking centre was 
set to a 20 oo =10 h 31 m 57.3 s and <5 2 ooo=-18°46'34.0" taken 
from the Nasa/IPAC Extragalactic Database (NED). The 
total integration time on source (6 antennae) amounted to 
^5 hours for the 2002 data and to ~3 hours for the 2003 
data. 

3. The data - results 

Both, the CO(1-0) and CO(2-l) lines are clearly detected 
(Fig.[2l[3]and[4]) while only upper limits can be determined 
for the continuum emission at 3 mm and 1 mm. 



M. Krips et al.: Barred CO emission in HE 1029-1831 



3 



Line 


line 


line 


integrated 


integrated intensity 


z 


M(H 2 +He) 6 




peak flux' 


width 


peak intensity 


in central —20" 








[mjy/beam] 


[km s- 1 ] 


[Jy/beamkm s _1 ] 


[Jy km s _1 ] 




[1O 9 X ] 


PdBI 


CO (1-0) 


65±7 


121±3 


8.4±0.2 


21±1 


0.0402±0.0001 


1-8 


CO (2-1) 


90±20 


131±7 


12±0.5 


25±1 


0.0402±0.0001 




BIMA 


CO(l-O) 


190±30 


150±20 


28±3 


30±3 


0.0402±0.0002 


2-12 




Centre of CO 


emission 


Dynamical Centre of CO emission 


PA C 


i d 




RA e 


Dec e 


RA e 


Dec 6 


[°] 


[°] 


PdBI 


CO (1-0) 


10 h 31 m 57.3 s - 


18°46'33.2" 


10 h 31 m 57.3 s 


-18°46'33.4" 


-20 


~ 10-30 


CO(2-l) 


10 h 31 m 57.3 s - 


18°46'33.2" 


10 h 31 m 57.3 s 


-18°46'33.4" 


-3 




BIMA 


CO(l-O) 


10 h 31 m 57.3 s - 


18°46'33.8" 


10 h 31 m 57.3 s 


-18°46'33.8" 


-0 


-10-30 



Table 1. Line parameters for CO(1-0) and CO(2-l) at the peak and integrated over the central 10-20". a taking 
a luminosity distance of 180 Mpc (based on the following cosmology: H = 70km s -1 , f2 m =0.3 and fix— 0.7; 1" 
corresponds then to — 0.8kpc). 6 assuming a M(Ha) to L' co conversion factor of 0.8-4.8 A4q (Kkm s _1 pc 2 ) -1 
(Solomon & Barrett 1991; Downes & Solomon 1998) and M(H 2 +Hc) = 1.36xM(H 2 ). c Estimated position angle of the 
CO emission (from North to East). d Estimated inclination angle of the CO disk. e The position uncertainties are of the 
order of -0.2" for the PdBI data and of -0.4" for the BIMA data. f flux errors include —10% (—20%) uncertainties 
from the flux calibration at 3 mm (1 mm). 
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Fig. 3. CO(1-0) (upper panel) and CO(2-l) (lower panel) 
line spectrum taken at the peak emission. The (optical) 
redshift of z=0.0403±0.0001 is indicated. The velocity res- 
olution is 3.4km s _1 (34 km s" 1 ) at both lines for the 
PdBI (BIMA). 



3.1. CO emission 

CO(1-0) emission was first detected with the BIMA ar- 
ray (Fig. [2]) and confirmed in subsequent observations 
of CO(1-0) and CO(2-l) with the IRAM PdBI. Figj3] 
shows the spectrum of the CO(1-0) and CO(2-l) line 
taken at the centroid of the CO emission for the IRAM 
PdBI and BIMA data. The redshift of the CO emission 
is determined to be z=0.0402±0.0001 for both transitions 
and both arrays which is consistent with the optical red- 
shift of z=0.0403±0.0001 (Fig. H Kaldare et al. 2003). 
A gaussian profile has been fitted to both lines. The cal- 
culated line parameters are given in Table [T] The inte- 
grated maps are plotted in Fig. [5] & HI Both CO lines are 
centered at the same position of aj2ooo =10:31:57.3 and 
Snooo = -18:46:33.2 (Table Q}. The CO(l-O) emission 
is not resolved by the modest angular resolution (<14") 
of the BIMA observations; the velocity integrated flux 
at peak and in the inner —20" are almost identical (see 
Tabled]). This excludes the existence of significant, very 
large scale (>10kpc) CO emission. The CO emission is, 
however, clearly extended in both transitions at the higher 
angular resolution (1-3") of the PdBI and reveals elonga- 
tions in north to south direction with a position angle 
(PA) of -0-20° (see Table [J) which is slightly smaller than 
the one derived for the optical bar (PA~30°; see Fischer 
et al. 2006). A combination of the BIMA and IRAM 
PdBI CO (1-0) maps enables an estimate of the CO size 
(=FWHM) by fitting an elliptical component to the uv 
data. We find a FWHM of -(7±2)" which corresponds to 
— (6±2) kpc. The integrated CO(1-0) line intensity results 
in a total gas mass of M(H 2 +He)«8 x 10 9 Ar© for the PdBI 
map and of — 1.2xl0 10 7V4© for the lower angular resolu- 
tion BIMA map assuming a M(H 2 ) to L' co conversion 
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Fig. 4. Upper panel: Velocity integrated CO(l-O) (left) 
and CO(2-l) [right) emission in HE 1029-1831 (be- 
tween ±200km s _1 ). Contours are from 3<r =0.8 Jy/beam 
km s _1 (1.4 Jy/beam km s _1 ) to 7.3 Jy/beam km s" 1 (13 
Jy/beam km s" 1 ) in steps of 3cr for CO(l-O) (CO(2-l)). 
Lower panel: Iso- velocity maps of the CO(l-O) and CO(2- 
1) emission CO(l-O) (left) and CO(2-l) (right) emission 
in HE 1029-1831. Contours are around the dynamic centre 
in steps of 10 km s _1 . The dashed lines indicate the cuts 
along which the position- velocity diagrams in Fig. [6] were 
taken. 



factor of 4.8 A4© (Kkm s 1 pc 2 ) 1 (Solomon & Barrett 
1991) and M(H 2 +He) = 1.36xM(H 2 ). The difference in in- 
tegrated intensity between the PdBI and BIMA maps in- 
dicates that the PdBI has filtered out approximately 30% 
of the total CO emission due to the lack of short-spacings. 
Taking the IRAS flux of 3.7 Jy at 3 THz (taken from NED) 
and assuming optically thin dust emission with a (rest- 
frame) dust temperatur^]] of ~40 K, we obtain a dust mass 
of ^lxlO 7 A^0. For a standard galactic gas-to-dust mass 
ratio of 150, the expected gas mass would be ~2xlO 9 A40, 
i.e., a factor of ^4-6 lower than using the CO luminosity. 
However, the assumed gas-to-dust mass ratio contains a 
large uncertainty and can increase in the cases of a few ac- 
tive Ultra-Luminous InfraRed Galaxies (ULIRGs) even up 
to 1000 or more (e.g., Contini & Contini 2003) so that the 
galactic value might not be applicable for HE 1029-1831. 
Alternatively, Downes & Solomon (1998) pointed out that 
the M(H2) to L' co conversion factor might be lower by a 
factor of 6 in the circumnuclear regions of ULIRGs. 



1 This value was estimated by fitting a grey body spectrum 
to the IRAS fluxes published by the NED. 
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Fig. 5. The CO(2-1)/CO(1-0) line ratio in temperature 
scale (color scale) overlaid with the velocity integrated 
CO(2-l) emission (contours). Increments for the velocity 
integrated CO(2-l) emission are the same as in Fig. 01 

As the angular resolutions of the PdBI for both 
transitions are comparable when using uniform weight- 
ing for CO(1-0) and natural weighting for CO(2- 
1), we can estimate the central (<8") line ra- 
tio i?2i=CO(2-l)/CO(l-0) (in temperature scale; i.e., 

R21 =(Sc02i(Jy)^ ob 2 S!C 02i)/( s coio)(Jy)^ "b s ,coio) for 
equal beamsizes) between both lines to be ~0.5, indicating 

subthermal excitation conditions and cold gas. However, 
one has to keep in mind that this estimate might still 
be biased by artifacts due to different uw-coverages for 
both lines in terms of spatial frequencies. One can reduce 
such a bias by apodising and truncating the respective 
itw-coverages to the overlapping region. This results in the 
map of line temperature ratios shown in Fig 5. The ratio 
map in this figure is not significantly different from one 
obtained simply from the unchanged uniformly weighted 
CO (1-0) and naturally weighted CO (2-1) maps. However, 
the resulting line ratio map still has to be interpreted with 
caution as it heavily depends on the adopted weighting 
factors. Thus, the integrated line ratio must be regarded 
as a best estimate. Although most of the CO appears to 
be subthermally excited and cold, the CO(2-l) line also 
indicates asymmetric emission at the south which has only 
a very weak counterpart in the CO(1-0) line. Such a dif- 
ference cannot be explained by resolution effects but could 
still be accounted for by a residual calibration error in the 
visibilities and a poor coverage of the uv-plane although 
no obvious inconsistencies were found. The high line ratio 
at the southern part of ~1 might thus indicate different 
excitation conditions of the gas. The location of the south- 
ern CO component is close to the southern spiral arm seen 
in the optical (see Fig. [I]). This might suggest that the CO 
at this position traces the contact point of the bar and the 
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southern spiral arm. The high line ratio could be then in- 
dicative of gas compression caused by a crossing of the 
bar- with the spiral arm orbits outside the bar corotation, 
consistent with standard theory. However, the southern 
CO emission has to be verified by further observations us- 
ing higher transitions of the CO line or denser gas tracers 
such as 13 CO or HCN. 

We have taken position velocity diagrams at all po- 
sition angles with a separation of 10° showing that the 
cuts in the direction of the CO extension (position an- 
gle: PA=0°) and perpendicular to it (PA=90°) represent 
the extremes in velocity structure: the cuts at PA=90° in 
both transitions reveal a clear and strong velocity gradi- 
ent (Fig. [H upper panels) which extends over a positional 
range of ±0.5"; by comparison, velocity gradients along 
the CO elongation, i.e., at PA=0°, are weaker and much 
less apparent (Fig. lower panels). Assuming a radius of 
~ 1" equivalent to 0.8 kpc and a velocity of ~ 80km s _1 , 
the total dynamical mass within the central 0.8 kpc - not 
corrected for inclination effects - can be calculated to be 
> 2 x lO 9 ^©. Probably, the gas disk is seen closer to 
face-on than to edge-on (compare also Fischer et al. 2006) 
so that this value must be regarded as lower limit. Since 
the dynamical mass includes the mass of the gas and of 
the stars, we can estimate the inclination of the galaxy 
by assuming that the gas mass of ~ 0.2 — 1.2 x 10 10 A4q 
contributes 10% to the total mass of the galaxy. This re- 
sults in an inclination of ~5-10°. By equating directly the 
estimated gas mass with the dynamical mass, the inclina- 
tion would increase to < 15-30° giving a more conservative 
upper limit. 
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Fig. 6. The upper panels show the position velocity di- 
agrams taken across the bar (PA=90°), and the lower 
panel those of a cut perpendicular to it, i.e., along the 
bar (PA=0°; see Fig. HJ. Increments are in steps of 10% 
from 20% and 30% of the maximum for the CO(1-0) and 
the CO(2-l) line respectively. The thick lines indicate the 
steepness of the gradient at PA=90° and are supposed to 
guide the reader's eyes. 



3.2. Continuum emission 

The continuum emission at 3 mm and 1 mm was not de- 
tected. The upper limits were determined in the line free 
channels (i.e., at velocities lower than —250km s _1 and 
larger than +250km s _1 at 3 mm and at velocities lower 
than —200km s" 1 and larger than —200km s _1 at 1 mm) 
to 3cr=3 mjy at 3 mm and to 3er = 8 mJy at 1 mm. The 
two upper limits are compatible with the mm-values pre- 
dicted by the grey body spectrum of the (infrared) dust 
emission (compare previous Section) which are of the or- 
der of ~l-5 mJy. 

4. Gas dynamics 

The kinematics and distribution of the CO emission in 
HE 1029-1831 suggest a barred potential in this galaxy 
(compare Roberts et al. 1979) as visible in the optical 
(Fischer et al. 2006). We have simulated the gas emis- 
sion in HE 1029-1831 by adopting the bar approach from 
Schinnerer, Eckart & Tacconi (2000) and Telesco & Decher 
(1988). Instead of computing the gas motions in a given 
gravitational potential, the gas is simply split up into sev- 
eral closed orbits with continous, non-overlapping curves 
of ellipticities e(r) and position angles (pa(r)) which are 
both smooth functions of the radius. Each ring, i.e., each 



radius has also an assigned velocity. For simplification, 
we adopt a constant rotation curve and a constant inten- 
sity over all rings. As fixed input parameters, we used the 
estimated PA of the emission of ^0-20° and the inclina- 
tion of ^10-30°. To find the best-fit model, we compared 
the simulated distribution and kinematics for different sets 
of e(r) and pa(r) with those from the observed gas emis- 
sion. Fig. [7] shows the result of these simulations. To fa- 
cilitate a comparison to the observed CO (1-0) data, the 
simulated CO distribution has been transformed to a uv- 
table based on the uv-coverage from the CO(1-0) PdBI 
observations; this is a more accurate method than just 
convolving the model with the beamsize of the CO(1-0) 
observations as it considers also spatial filtering effects 
due to the discrete uv-coverage. This simplified model al- 
ready reproduces the observed emission. The simulated 
distribution and the kinematics of the gas agree very well 
with each other indicating that the gas properties can be 
almost completely explained by assuming a barred poten- 
tial. However, the asymmetry between the CO(1-0) and 
CO(2-l) line cannot be explained by the symmetric bar 
model. 

2 The used code 3DRings is described in more detail by 
Schinnerer, Eckart & Tacconi (2000). 
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Fig. 7. Upper panel: Iso- velocity diagrams of the CO(l- 
0) emission in HE 1029-1831 (left) and of a first-order 
bar model (right) overlaid to the respective integrated 
line emission maps. Increments are again in steps of 
10 km s _1 .Lower panel: Position-velocity diagram simu- 
lated with a bar model for the emission in HE 1029-1831 
taken along a cut with a PA of 90° (left) and of 0° (right); 
compare Fig. 2] & [SI The dashed contours are superim- 
posed on the observed data presented in Fig.[6j Increments 
are in steps of 10% from 20% of the maximum. 



5. Summary and Conclusions 

Strong and extended CO(1-0) and CO(2-l) emission is 
detected in the barred Seyfert galaxy HE 1029-1831 with 
the IRAM PdBI and the BIMA array while no contin- 
uum emission was found. Both lines agree well with each 
other in distribution, position and velocity width and im- 
ply a total gas mass of ~ 1.2 x lO^A'f©. The difference 
between the dynamical mass (~2xlO 9 A^0, not corrected 
for inclination) indicates a low inclination of ^10-30°. The 
CO(1-0) and CO(2-l) emission are extended along the op- 
tical bar. The FWHM of the CO emission is calculated to 
be ~(6±2)kpc. Also, a striking velocity gradient is found 
at a position angle of 90°, i.e., along the minor axis of 
the bar, extending over ~ 1" approximately. In the direc- 
tion of the bar (PA=0°), almost no velocity gradient is 
seen. Adopting the approach from Schinnerer, Eckart & 
Tacconi (2000), the line emission can be simulated with 
a very simple bar model explaining the detected distribu- 
tion and kinematics. Following Roberts et al. (1979), the 
strong velocity gradient across the bar (PA=90°) can be 
interpreted as bar-driven inflow of the gas while the mo- 
tions along the bar (PA=0°) indicate a velocity dispersion 



of Av=±80km s -1 . The line ratio between CO(2-l) and 
CO(1-0) is estimated to be ^0.5 over most areas of the 
bar although probably still biased by resolution effects. 
This indicates subthcrmally excited and cold gas which is 
typical for gas in a bar/disk. The southern increase of the 
CO(2-1)/CO(1-0) line ratio to over 1 indicates different 
excitation conditions which might be caused by gas com- 
pression at the crossing point of the bar and the southern 
spiral arm. However, this southern CO component still 
requires further support. 
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